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The BABAR detector, which operates at the SLAC PEP-II asymmetric e + e~ collider 
at energies near the T(4S) resonance, started operation on the 26 th of May 1999. We 
present the first study of sin2/5, with samples of B° — > J/ipK® and B° — > ip{2S)K® decays, 
using 9.0 ftT 1 of data recorded between January and July 2000 at the Y(AS) resonance 
and 0.8 fb" 1 recorded 40MeV below the T(4S) resonance. A preliminary result of sin2/3 = 
0.12 ± 0.37 (stat) ± 0.09 (syst) was obtained. Details of the analysis are given. Moreover, 
we present measurements of charged and neutral B meson lifetimes and B°B° oscillation 
frequency. 
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1 Introduction 



The primary goal of the BABAR experiment at PEP-II is to overconstrain the Unitarity 
Triangle. The sides of this triangle can be measured through non-CP violating physics, 
such as V u b, V c b, Vtd measurements [0], while its angles are accessible through CP violating 
processes [Q. 

2 PEP-II 

The PEP-II B Factory is an e + e~ colliding beam storage ring complex on the SLAC site 
designed to produce a luminosity of at least 3xl0 33 cm _2 s~ 1 at a center-of-mass energy of 
10.58 GeV, the mass of the T(4S) resonance. In the 2000 run, the achieved average luminosity 
was 2.5xl0 33 cm _2 s _1 , with a daily average integrated luminosity of 135 pb - . The total 
collected luminosity was about 22 fb~ . The machine is asymmetric with a High Energy Ring 
(HER) for the 9.0 GeV electron beam and a Low Energy Ring (LER) for the 3.1 GeV positron 
beam. This corresponds to ,$7=0. 56 and makes it possible to measure time dependent CP 
violating asymmetries. It corresponds to an average separation of /57cr=250/im between 
the two B mesons vertices. 

3 BaBar 

3.1 Detector description |2|] 

The volume within the BABAR superconducting solenoid, which produces a 1.5 T axial mag- 
netic field, consists of: a five layer silicon strip vertex detector (SVT), a central drift chamber 
(DCH), a quartz-bar Cherenkov radiation detector (DIRC) and a Csl crystal electromagnetic 
calorimeter (EMC). Two layers of cylindrical resistive plate counters (RPCs) are located be- 
tween the barrel calorimeter and the magnet cryostat. All the detectors located inside the 
magnet have full acceptance in azimuth. The integrated flux return (IFR) outside the cryo- 
stat is composed of 18 layers of steel, which successively increase in thickness away from the 
interaction point, and are instrumented with 19 layers of planar RPCs in the barrel and 18 
in the endcaps. 

3.2 Event reconstruction j2| 

Charged particles are detected and their momentum is measured by a combination of 
the DCH and SVT. The charged particle momentum resolution is approximately given by 
(Spt/pt) = (0.0015j9 T ) 2 + (0.005) 2 , where Pt is in GeV/c. The SVT, with a typical res- 
olution of 10 /im per hit, provides excellent vertex resolution both in the transverse plane 
and in z. The vertex resolution in z is typically 50 /im for a fully reconstructed B meson 
and of order 100 /im for the distance among the two B mesons when only one is fully re- 
constructed. Leptons and hadrons are identified using a combination of measurements from 
all the BABAR components, including the energy loss dE/dx in the helium-based gas of the 
DCH (40 samples maximum) and in the silicon of the SVT (5 samples maximum). Electrons 
and photons are identified in the barrel and the forward regions by the EMC, and muons are 



identified in the IFR. In the barrel region the DIRC provides excellent kaon identification 
over the full momentum range above 250MeV/c. 



4 sin2/3 measurement 

In e + e~ storage rings operating at the T(45) resonance a B°B° pair produced in a T(45) 
decay evolves in a coherent P-wave until one of the B mesons decays. If one of the B 
mesons (B tag ) can be ascertained to decay to a state of known flavor at a certain time t tag , 
the other B (Bcp) is at that time known to be of the opposite flavor. For the measurement 
of sin2/3, Bqp is fully reconstructed in a CP eigenstate (J/ipK® or ip(2S)Kg). By measuring 
the proper time interval At = tcp — ttag fro m the B tag decay time to the decay of the Bqp 
(tcp), it is possible to determine the time evolution of the initially pure B° or B° state: 

/±( At; T, Am d , £> sin 2/3) = ire~ r|At| [1 ± £> sin 2/3 x sinAm d At] , (1) 

where the + or — sign indicates whether the B tag is tagged as a B° or a B°, respectively. 
The dilution factor T> is given by D = 1 — 2w, where w is the mistag fraction, i.e., the prob- 
ability that the flavor of the tagging B is identified incorrectly. A direct CP violation term 
proportional to cos Am,i At could arise from the interference between two decay mechanisms 
with different weak phases. In the Standard Model, we consider that the dominant diagrams 
for the decay modes have no relative weak phase, so no such term is expected. 

To account for the finite resolution of the detector, the time-dependent distributions f± 
for B° and B° tagged events (Eq. [I]) must be convoluted with a time resolution function 
TZ(At;d): 

F±( At; T, Am d , £>sin2/3,a) = /±(At; T, Am d , V sin2/3 ) ® K{ At ; a) , (2) 

where a represents the set of parameters that describe the resolution function. 
It is possible to construct a CP-violating observable 

T+{At) — -F-(At) 
Acp{M) ~ f + (At) + rjAl) ' (3) 

which is approximately proportional to sin2/?: 

A C p{At) ~ V sin 2/3 x sin Am d At . (4) 

Since no time-integrated CP asymmetry effect is expected, an analysis of the time- dependent 
asymmetry is necessary. At an asymmetric-energy B Factory, the proper decay-time differ- 
ence At is, to an excellent approximation, proportional to the distance Az between the two 
P°-decay vertices along the axis of the boost, At rs Az/c (/3j). 

Since the amplitude of the time- dependent CP-violating asymmetry in Eq. [| involves 
the product of V and sin2/3, one needs to determine the dilution factors T>i (or equivalently 
the mistag fractions Wi) in order to extract the value of sin2/3. The mistag fractions are 
determined from the data by studying the time-dependent rate of B°B° oscillations. 




Beam-energy substituted mass (GeV/c ) 



Beam-energy substituted mass (GeV/c ) 



Figure 1: J/ipK° s (K® -> tt+tt-J si^na/ (left). J/ipK° s (K° s -> vrVj si 5 na/ fr^faj. 



The value of the single free parameter sin2/3 is extracted from the tagged Bqp sample by 
maximizing the likelihood function 



CP 



E 



J2 ln F+{ At; T, Am d) a, sin 2/3) + ^ lnjF_( At; T, Am d , a, sin 2/3) 

B° tag £0 tag 

(5) 

where the outer summation is over tagging categories i. 



4.1 Analysis 

For this analysis we use a sample of 9.8 fb" 1 of data recorded between January and the 
beginning of July 2000, of which 0.8 fb -1 was recorded 40MeV below the T(4S) resonance 
(off-resonance data). 

The measurement of the CP-violating asymmetry has five main components : 

• Selection of the signal B°/B° _> J/ijjK 3 and B°/B° -> ip(2S)K° events, as described 
in detail in 0. 

Distributions of AE, the difference between the reconstructed and expected B meson 
energy measured in the center-of-mass frame, and mss, the beam-energy substitued 
mass, are shown in Fig. [I] for the J/ipK° s (K° s — > n + ix~) and JfyK® (K® — > 7r°7r°) 
samples. Signal event yields and purities, determined from a fit to the uies distributions 
after selection on AE, are summarized in Table ID. 



Measurement of the distance Az between the two B° decay vertices along the T(4S) 
boost described in detail in j|] and 0. 



Table 1: Event yields for the different samples used in this analysis, from the fit to toes distribu- 
tions after selection on AE. The purity is quoted for toes > 5.270 MeV/c 2 . 



Final state 


Yield 


Purity (%) 


Jft>K% (K [ > - 


-> 7T + 7T ) 


124±12 


96 


J/^K° S (K° s - 


- vr vr ) 


18±4 


91 






27±6 


93 



From the measurement of Az , At can be computed. 

The time resolution function is described accurately by the sum of two Gaussian dis- 
tributions, which has five independent parameters: 

K(At;a) = £ — ^= exp (-(At - <y,)7W) . (6) 

A fit to the time resolution function in Monte Carlo simulated events indicates that 
most of the events {j\ = 1 — ji = 70%) are in the core Gaussian, which has a width 
0i ~ 0.6 ps. The wide Gaussian has a width cr 2 rs 1.8 ps. Tracks from forward-going 
charm decays included in the reconstruction of the B tag vertex introduce a small bias, 
Si ?z — 0.2ps, for the core Gaussian. 

• Determination of the flavor of the B tag , as described in detail in 0. 

Each event with a CP candidate is assigned a -B or B° tag if the rest of the event (i.e., 
with the daughter tracks of the Bcp removed) satisfies the criteria from one of several 
tagging categories. In other words, a B° tag indicates that the Bcp candidate was in 
a B° state at At = 0; a B° tag indicates that the Bqp candidate was in a B° state. 

Three tagging categories rely on the presence of a fast lepton (Lepton categories) 
and/or one or more charged kaons in the event (Kaon category). Two categories, called 
neural network categories (NT1 and NT2), are based upon the output value of a neural 
network algorithm applied to events that have not already been assigned to lepton or 
kaon tagging categories. 

• Measurement of the dilution factors T>i from the data for the different tagging cate- 
gories, as described in detail in j4|. 

The figure of merit for each tagging category is the effective tagging efficiency Qi = 
Si (1 — 2wi) 2 , where is the fraction of events assigned to category % and Wi is the 
mistag fraction. 

The mistag fractions are measured directly in events in which one B° candidate, called 
the B rec , is fully reconstructed in a flavor eigenstate mode. The flavor-tagging al- 
gorithms are applied to the rest of the event, which constitutes the potential B tag , 
assuming the B rec is properly tagged. 

The mistag fractions and the tagging efficiencies obtained by combining the results from 
maximum likelihood fits to the time distributions in the B° hadronic and semileptonic 
samples are summarized in Table 0. 



Table 2: Mistag fractions measured from a maximum-likelihood fit to the time distribution for the 
fully-reconstructed B° sample. The Electron and Muon categories are grouped into one Lepton 
category. The uncertainties on e and Q are statistical only. 



Tagging Category 


£(%) 


w (%) 


Q{%) 


Lepton 


11.2 ±0.5 


9.6 ± 1.7 ± 1.3 


7.3 ±0.7 


Kaon 


36.7 ±0.9 


19.7 ± 1.3 ± 1.1 


13.5 ±1.2 


NT1 


11.7 ±0.5 


16.7 ±2.2 ±2.0 


5.2 ±0.7 


NT2 


16.6 ±0.6 


33.1 ±2.1 ±2.1 


1.9 ±0.5 


all 


76.7 ±0.5 




27.9 ± 1.6 



Table 3: Summary of systematic uncertainties. The different contributions to the systematic error 
are added in quadrature. 



Source of uncertainty 


Uncertainty on sin2/3 


uncertainty on Tg 


0.002 


uncertainty on Am^ 


0.015 


uncertainty on Az resolution for CP sample 


0.019 


uncertainty on time-resolution bias for CP sample 


0.047 


uncertainty on measurement of mistag fractions 


0.053 


different mistag fractions for CP and non-CP samples 


0.050 


different mistag fractions for B° and B° 


0.005 


background in CP sample 


0.015 


total systematic error 


0.091 



• Extraction of the amplitude of the CP asymmetry and the value of sin2/3 with an 
unbinned maximum likelihood fit. 

The maximum-likelihood fit for sin2/3, using the full tagged sample of B°/B° — ► J/if)K° s 
and B°/B° -> ^(2S)K° S events, gives: 

sin2/3 - 0.12 ± 0.37 (stat) ± 0.09 (syst). (7) 

For this result, the B° lifetime and Arrid are fixed to the current best values 0. The log 
likelihood is shown as a function of sin2/3 and the At distributions for B° and B° tags 
are shown in Fig. ^. The contributions to the systematic uncertainty are summarized 
in Table |3]. They are added in quadrature to obtain the total systematic error. 

Improvements in the sin2/5 measurement are expected in the near future with the 
accumulation and analysis of more data and further systematic studies. 
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Figure 2: Variation of the log likelihood as a function of sin2/3 (left). The two horizont al dashed 
lines indicate changes in the log likelihood corresponding to one and two statistical standard devi- 
ations. Distribution of At for (a) the B° tagged events and (b) the B° tagged events in the CP 
sample (right). 



5 Measurements of charged and neutral B meson lifetimes and 
.B ^ oscillations 

These measurements can be used to test theoretical models of heavy-quark decays and to 
constrain the Unitarity Triangle (via the sensitivity to the value of the CKM matrix element 

v td ). 

One B (B rec ) is fully reconstructed in an all-hadronic (B° — > D ( ~*^n + , p + , D^^af, 

J/ipK*° and B + -> D { * ] \ + , J/ipK+, ip(2S)K+) or semileptonic decay mode (B° -> D*-£+u) Q. 
A total of about 2600 B° and a similar number of charged B are reconstructed in the hadronic 
modes, with a mean purity of ~ 90%. The background is mainly combinatorial. About 7500 
B° are reconstructed in the semileptonic modes, with a mean purity of ~ 84%. Backgrounds 
to the semi-leptonic mode are due to combinatorics, wrong leptons, cc events, and charged 
B decays from B + — > D**°l + u. 

The separation between the two B vertices (z rec and z?ag f° r t ne reconstructed and 
tagged vertex, respectively) along the boost direction, Az = z rec — z tag , is measured and 
used to estimate the decay time difference, At = Az/(3 z c. The At resolution is dominated 
by the precision on the B tag vertex, and has little dependence on the decay mode of the B rec . 



1 Throughout this paper, conjugate modes are implied. 
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Figure 3: At distributions for B°/B° (left) and B + /B~ (right) candidates in the signal region 
(mES> 5.27 GeV/c 2 ). The result of the lifetime fit is superimposed. The background is shown by 
the hatched area. 

5.1 Lifetime Measurements 

The B° and B + lifetimes are extracted from a simultaneous unbinned maximum likelihood 
fit to the At distributions of the signal candidates, assuming a common resolution function. 
Only hadronic modes have been used. An empirical description of the At background shape is 
assumed, using m^s sidebands with independent parameters for neutral and charged mesons. 
Fig. |H shows the At distributions with the fit result superimposed. 

5.2 Time— dependent B°B° mixing 

A time-dependent B°B° mixing measurement requires the determination of the flavor of both 
B mesons. Considering the B°B° system as a whole, one can classify the tagged events as 
mixed or unmixed depending on whether the B tag is tagged with the same flavor as the B rec 
or with the opposite flavor. 

From the time-dependent rate of mixed (N m i X ) and unmixed (N unm i X ) events, the mixing 
asymmetry a(At) = (N unmix - N mix ) / \N unmix + N mix ) is calculated as a function of At 
and fit to the expected cosine distribution. A simultaneous unbinned likelihood fit to the 
mixing asymmetry frequency and its amplitude allows the determination of both Am^ and 
the mistag rates, Wi, for the different tagging categories. The fit is performed simultaneously 
in each tagging category, assuming a common resolution function. Fig. f| shows the a(At) 
distributions with the fit result superimposed. 
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Figure 4: Time- dependent asymmetry a(At) between unmixed and mixed events for (left) hadronic 
B candidates with mgs > 5.27 GeV/c 2 and (right) for B — > D*lv candidates. 

5.3 Results 

The preliminary results for the B meson lifetimes are t b o = 1.506±0.052 (stat)±0.029 (syst) ps 
and t b + = 1.602 ± 0.049 (stat) ± 0.035 (syst) ps and for their ratio is t b +/t b o = 1.065 ± 
0.044 (stat) ± 0.021 (syst). 

From the hadronic B° sample we measure the B°B° oscillation frequency: Amd = 0.516 ± 
0.031 (stat) ± 0.018(syst) ftps -1 and from the D*~£ + u sample the result is Am d = 0.508 ± 
0.020 (stat) ± 0.022(syst) fops -1 Combining the two Am^ results, we obtain the preliminary 
result: Am d = 0.512 ± 0.017(stat) ± 0.022(syst) hps -1 . The mistag rates and At resolution 
function extracted from these fits to the data are used in the BABAR CP violation asymmetry 
analysis 0]. The above results are consistent with previous measurements and are of 
similar precision. The mixing result is compatible with a BABAR measurement using di- 
leptons 0. 

6 Conclusions 

We have presented BABAR' s first measurement of the CP- violating asymmetry parameter 
sin2/? in the B meson system: 

sin2/? = 0.12 ± 0.37 (stat) ± 0.09 (syst). (8) 

Our measurement is consistent with the world average sin2/3 = 0.9 ±0.4 ||, and is currently 
limited by the size of the CP sample. 

We have presented also the time-dependent mixing and lifetime measurements, performed 
for the first time at the T(AS). 



Other important competitive measurements have been achieved by BABAR. 
All measurement are foreseen to be improved in the near future, with the accumulation 
and study of more data and further systematic studies. 
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